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JOULE -THOM SON INVERSION CURVES AND RELATED 
COEFFICIENTS FOR SEVERAL SIMPLE FLUIDS 
by Robert C. Hendricks, Ildiko C. Peller, and Anne K. Baron 

Lewis Research Center 

SUMMARY 

The pressure- volume-temperature (PVT) relations for several fluids (methane, 
nitrogen, oxygen, argon, carbon dioxide, neon, carbon monoxide, fluorine, hydrogen, 
and helium) have been used to study the behavior of the isenthalpic Joule- Thomson coef- 
ficient and in particular the inversion curve. Some attention is also given to Joule- 
Thomson refrigerators and liquifiers. For the fluids methane, nitrogen, oxygen, argon, 
carbon dioxide, and neon the calculations involved extrapolation of the equation of state 
to high pressure and temperature. In all cases the PVT extrapolations and the deriva- 
tives (SP/aT)^ and (3P/3 p),p, where p is density, appeared to be natural extensions of 
the lower pressure and temperature results. These results for nitrogen were com- 
pared to the PVT relation of Coleman and Stewart, which was developed for nitrogen and 
covers pressures to 1000 atmospheres and temperatures to 1000 K. The comparisons 
were quite favorable. 

From these results a Joule- Thomson inversion curve was calculated for each fluid. 
The curves were compared to existing data with favorable results for all but carbon 
monoxide, where sizable deviations from the data were noted for temperatures above a 
reduced temperature T^ of 3. An approximation for the critical isenthalpic Joule- 
Thomson coefficient p c was developed. This simple formula, p c = T c /6P c , was found 
to be in good agreement with the data of Roebuck and Osterberg and Coleman and Stewart 
for nitrogen. Reduced Joule- Thomson data were determined for eight fluids and plotted 
as a function of reduced temperature for various isobars (reduced pressures). 

The inversion curve and the isothermal Joule-Thomson coefficient are shown to 
comply with the principle of corresponding states. The isenthalpic Joule-Thomson var- 
iable does not. Nevertheless, the reduced inversion curves are useful when data are 
absent and for grouping data in fluid- to- fluid comparisons. 

The inversion temperature and the Boyle temperature as calculated from the second 
virial coefficient B were also found to be in agreement with values found in the litera- 
ture. 

Finally, a revised form of the generalized inversion curve is presented. The fluids 
helium, hydrogen, and neon were not included in the generalization because of their 
quantum effects. Carbon dioxide was also excluded from the generalization because of 
its high acentric factor. Available low- temperature- fluorine data were included. 



INTRODUCTION 


Joule- Thomson expansion is widely used for liquification of gases and for refrigera- 
tion. There are no moving parts, maintenance is minimal, and reliability is high. 

Joule- Thomson expansion is usually used in the lowest temperature portion of a liquifi- 
cation or refrigeration cycle; however, when coupled with a pump, a complete cycle is 
formed. The Joule- Thomson liquification and refrigeration processes are discussed in 
appendix D along with some indication of how the figure of merit (cost of producing a 
unit of liquid or cooling) influences the cycle design. It should be noted that the maxi- 
mum refrigeration or liquification does not correspond to minimum cost. A variety of 
specific cooling and liquification applications can be found in reference 1. 

Under certain conditions a Joule- Thomson expansion will produce heating rather 
than cooling; consequently, determining these conditions is important to the design of 
cryogenic equipment. In the Joule- Thomson experiment, the gas is expanded isenthal- 
pically through a porous plug, and the slope of the isenthalpic curve (constant H) in the 
pressure- temperature (P-T) plane is called the Joule- Thomson coefficient /i 



( All symbols are defined in appendix F. ) The unique locus of points where the isenthalp 
is a maximum is called the inversion curve. Thus, for a certain temperature, there 
exists a pressure beyond which ju < 0 and isenthalpic expansion causes a temperature 
rise. However, at lower pressures, fi > 0 and cooling results from isenthalpic expan- 
sion. The maximum isenthalpic cooling occurs when conditions on the high-pressure 
side of the plug (or valve) fall on the inversion locus (appendix D). 

While the Joule- Thomson effect has been experimentally determined for several 
fluids (e.g. , refs. 2 to 7), it is not only a very difficult experiment but requires precise 
PVT data for cross interpolation. Thus, when accurate experimental data for a parti- 
cular gas are either inadequate or nonexistent, it becomes necessary to estimate the 
inversion curve. While the principle of corresponding states has been applied to the 
inversion curve by several authors (e.g. , refs. 8 and 9), Gunn, Chueh, and Prausnitz 
(ref. 10) reexamined both the experimental and calculated Joule- Thomson "data" for 
argon, xenon, nitrogen, carbon monoxide, methane, ethylene, neon, helium, hydrogen, 
and some mixtures (see also appendix A). Quantum gases (helium, hydrogen, and neon) 
do not follow the corresponding states principle except at high temperatures and, con- 
sequently, were correlated separately. The inversion curve for the remaining gases 
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Figure 1. - Generalized reduced inversion curve for several corre- 
sponding states fluids as determined by Gunn, Ctiueh, Prausnitz 
(ref. 10). 


(fig. 1) was well represented in terms of reduced temperature (Tj^ = T/T ) and re- 
duced pressure (Pp^ = P/P ) by 
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Gunn, Chueh, and Prausnitz (ref. 10) also concluded that Joule- Thomson coeffi- 
cients calculated with an accurate equation of state should be more accurate than those 
measured experimentally. It must be remembered that accuracy in measuring the pres- 
sure drop in an isenthalpic experiment and interpolating these data to determine 
(3T/3 P)jj is very difficult to obtain. 

Recently, accurate PVT data on methane and oxygen have been acquired (refs. 11 
and 12). Coleman and Stewart (ref. 13) have reexamined all the available nitrogen data 
and presented an equation of state that is valid over a considerable range of pressure. 
Bender (ref. 14) reported an accurate equation of state for five fluids: methane, nitro- 
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gen, oxygen, argon, and carbon dioxide. Hendricks, Baron, Peller, and Pew (ref. 15) 
developed a computer program to calculate the thermodynamic and transport properties 
for eight fluids. The primary equation selected in this program was the Bender equa- 
tion, because it covered five fluids of considerable interest with a consistent curve fit. 
Carbon monoxide, neon, and helium were also included in reference 15 by using equa- 
tions formulated at NBS- Boulder (see ref. 15). In addition, the Coleman and Stewart 
equation for nitrogen (ref. 13), because of its extended pressure range, was also in- 
cluded in the program of reference 15. 

For most fluids, neither the equation of state nor the data cover the extreme 
pressure-temperature ranges required to obtain the Joule- Thomson curve and related 
results. The only equation used in reference 15 which was originally fit to the very high 
pressures required to determine the inversion curve was the Coleman and Stewart equa- 
tion for nitrogen (ref. 13). Consequently, the first topic of this report is a discussion 
of the validity of the PVT results obtained when the equations of state for the eight fluids 
as defined in reference 15 are extrapolated to the high pressure required to compute the 
Joule- Thomson inversion curve. 

Second, by using the computer program of reference 15, the work of reference 10 
is extended to present inversion curves and related Joule- Thomson results for the fluids 
methane, nitrogen, oxygen, carbon dioxide, argon, carbon monoxide, neon, and helium. 
Limited results for hydrogen, fluorine, and some noble gases are also included. 

The next topic covered is the calculation of inversion temperatures by using the 
second virial coefficient of the equation of state programmed in reference 15. 

Finally, the principle of corresponding states as applicable to the inversion curve 
and the critical isenthalpic Joule- Thomson coefficient is discussed. 


HIGH-PRESSURE PVT EXTENS IONS 

In order to extrapolate Bender's equation of state (ref. 14) to higher pressures and 
temperatures, it is first necessary to examine the extended isochores. The numerical 
values of the extended isochores for each fluid of figure 2 were based on the following 
increments of the critical isochore: 0.5 p , 1.0 p., 1.25 p„, 1.5p„, 2.0p„, 2.5 p„, 
and 3. 0 p c . The isochores of figure 2 are bounded above by a pressure of 100 MN/ni 
and from the right by a temperature around 600 K. While these bounds may at first 
seem arbitrary, they were established to encompass most of the Joule- Thomson loci and 
in the case of nitrogen to compare the extrapolated isochores as predicted by Bender 
(ref. 14) to those predicted by Coleman and Stewart (ref. 13). 

The equations of state for the fluids methane, nitrogen, oxygen, argon, carbon diox- 
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1 2 
ide, carbon monoxide, and neon were all fitted for pressures less than 40 MN/m , and 

p 

that for helium for pressures less than 10 MN/m (ref. 15). Nevertheless, the extrapo- 
lated isochores appear linear, continuous, and a natural extension from the lower pres- 
sure regime. 

The slopes (3P/ST)^ as presented in figure 3 indicate the linearity of the isochores. 

In general, the (9P/9T) curves are nearly constant for densities to p/p <=* 2. How- 

P ^ 

ever, at higher densities, that is, p/p of 2.5 and 3.0, the (9P/9T) exhibits a high 

c p 

degree of curvature at low temperatures and tends toward a constant at higher temper- 
atures. 

The derivatives (3P/3p)rp, for various fluids, are not as easily compared. However, 
similar trends were noted for densities to p/p c of approximately 2 (perhaps 2.5), even 
though quantitative trends were not established. 

Due to the linear nature of the isochores in the P-T plane, the behavior of these 
derivatives is not unusual; and, as such, one would expect the extrapolated PVT results, 
in the absence of data, to be fairly good. A second reason for investigating the behavior 
of these derivatives is that they are used in the calculation of the Joule- Thomson coef- 
ficient. 

2 

The Coleman- Stewart equation of state for nitrogen (ref. 13) is valid to 100 MN/m 
and 1000 K (this region includes the inversion curve). The extrapolated values of Ben- 
der's equation in figures 2(b) and 3(b) for nitrogen should agree with those calculated by 
using the Coleman- Stewart equation, in order to have greater confidence in the extrapo- 
lated results. 

Apart from a slight shift of the high- density isochore as shown in figure 4, the iso- 
chores calculated by the two techniques are the same. In general, (9P/3T)^, as calcu- 
lated from references 13 and 14 by reference 15, are in good agreement (fig. 5). 

In view of these results, one can expect the extrapolated PVT results to be a reason- 
able approximation to high-pressure data (when it becomes available). 


^"Preliminary Joule- Thomson points for fluorine are also presented; however, due 
to some unexpected curvature in the isochores and deviations in the enthalpy calculation, 
the J-T locus was not calculated herein. 
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Figure 3. - Partial derivative (dP/dTL as calculated from PVT extensions for selected isochores (ref. 15), for various critical densities p c . 
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Figure 4. - Comparison of pressure-temperature plane for selected isochores 
based on reference 13 and the extended PVT results based on reference 14. 
Critical density, p c , 0. 3105 g/cm 3 . 
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Figure 5. - Comparison of partial derivative (dP/dT)p based on reference 13 and 
the extended PVT results based on reference 14. Critical density, p c , 0. 3105 
g/cm 3 . 



JOULE-THOMSON EFFECTS 


Inversion Curves 

The isenthalpic Joule- Thom son coefficient, defined as the slope of an isenthalp in 
the T-P plane, can be written as 


M = 




1 


(1) 


The inversion curve is defined as the locus of maxima of these isenthalps in the T-P 
plane. This locus of points, where ju = 0, can be calculated by determining the maxima 
or more simply from 



( 2 ) 


A discussion of the relation between the isenthalpic and isothermal Joule- Thomson 
coefficients and their relation to the principle of corresponding states is found in appen- 
dix A. As pointed out in appendix A, the isenthalpic Joule-Thomson coefficient jit is 
not a corresponding states variable. However, the isothermal Joule-Thomson coeffi- 
cient cp is a corresponding states variable. They are related by 

<P =-MC p (A3) 

Now when cp - 0 ju = 0 since Cp =£ 0; and consequently, the principle can be applied to 
the locus of points where ju = 0. 

The following discussion of inversion curves includes the nonreduced and reduced 
forms. Two fluids, nitrogen and oxygen, because they are well represented with re- 
liable PVT data and J-T data in the literature, are treated in both forms. Reduced 
curves are presented for the eight fluids of reference 15 along with some preliminary 
results for hydrogen and some noble gases. 
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Non reduced Inversion Curves for Nitrogen and Oxygen 


Nitrogen is a very useful research fluid and thus several sets of data for the inver- 
sion curve appear in the literature. The experimental Joule- Thomson inversion curve 
from the data of Roebuck and Osterberg (ref. 5), Deming and Shupe (refs. 16 and 17), 
Limbeck, Michels, and Wolkers (ref. 18), Dean and Mann (ref. 19), and Din (ref. 20), 
along with the inversion curve calculated from the equation of reference 10, are plotted 
in figures 6(a) and (b). The inversion curve as calculated by using the Coleman- Stewart 
equation of state (ref. 13) is illustrated as a set of data points in figure 6(a) in order to 
distinguish it from the inversion curve calculated by using the extrapolated Bender equa- 
tion of state (ref. 14), which is plotted as a dashed line. As can be seen, the agreement 
between the data and the calculated inversion curve is quite good, indicating that such 
calculations for other fluids should, in the absence of data, be reasonably good approxi- 
mations. Also plotted are selected enthalpies to illustrate their general behavior in the 
P-T plane and in particular near the inversion and saturation curves. 

In figure 6(b) the exercise of figure 6(a) is repeated by using the Coleman -Stew art 
equation of state (ref. 13) to generate the dashed line. As can be seen, the inversion 
curves are in good agreement, as are the selected isenthalps. This lends a great deal 
of confidence to the extrapolated PVT results. 

Oxygen has significant applications in the nation's space program and in industrial 
usage. The recent PVT data of Weber (ref. 12) are very accurate, and a portion of the 
inversion curve was defined by Weber (ref. 12). The inversion curve calculated by us- 
ing the extrapolated PVT results of Bender (ref. 14) and the data of Weber (ref. 12) are 
given in figure 7. As can be seen, the agreement is quite good; however, there is no 
verification of the high- temperature portion of the inversion curve. Again selected en- 
thalpy values are overlaid on figure 7 to illustrate isenthalpic behavior for oxygen in the 
P-T plane. 
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Reduced Inversion Curves 

In references 8 to 10, to cite a few, the principle of corresponding states was used 
to group the inversion curves. (See appendix A on the applicability of the principle. ) 
Furthermore, Gunn, Chueh, and Prausnitz (ref. 10) found that several fluids (nitrogen, 
argon, xenon, carbon monoxide, methane, and ethylene) could be represented by a 
single inversion curve, while quantum fluids such as neon, hydrogen, and helium would 
not fit the pattern. 

Since the principle of corresponding states should be applicable at high tempera- 
tures, Gunn, Chueh, and Prausnitz (ref. 10) force- fit the inversion data for hydrogen 
and heli um to make their maximum inversion temperature coincide with that of the other 
fluids, and thereby defined pseudocritical values for the two fluids (refs. 10 and 21). 

The inversion curves for methane, nitrogen, oxygen, argon, carbon dioxide, car- 
bon monoxide, neon, nitrogen (Coleman and Stewart, ref. 13), and helium were calcu- 
lated by using the extrapolated PVT results and plotted in reduced coordinates as fig- 
ures 8 and 9. The inversion curve for hydrogen (fig. 9(c)) and portions of the inversion 
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(a) Methane. P c ■ 4.627 MN/m 2 ; T c > 190.77 K. 


(b) Nitrogen. P c ■ 3.417 MN/m 2 ; T r - 126. 3 K. 




(c) Nitrogen (based on ref. 13; not a high pressure-temperature ex- 
tension) P c - 3.4 MN/m 2 ; T c " 126. 2 K. 


(dl Oxygen. P c • 5. 083 MN/m 2 ; T - 154. 78 K. 


Figure 8. - Reduced Joule-Thomson inversion curve for several fluids, at various critical pressures P c and critical temperatures T c based 
on PVT extensions as calculated from reference 15. Available data included for comparison. 
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Figure 9. - Reduced Joule-Thomson inversion curves for the quantum fluids helium, hydrogen and 
neon, for various values of critical pressure P c and critical temperature T c . 
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Figure 9. - Continued. 



curves for xenon and krypton are presented (figs. 10 and 11). The helium, hydrogen, 

and neon pseudoreduced plots as suggested by Gunn, Chueh, and Prausnitz (ref. 10) 

are also presented for comparison (figs. 9(b), (d), and (f)). The maximum inversion 

temperature (above which expansion produces heating at all pressures) as calculated in 

appendix B and summarized in table I is plotted on the inversion curve for each fluid. 

The Boyle temperature T B and the low-density Joule inversion temperature T B 

a max 

are also discussed in appendix B and tabulated in table I. However, for fluids consid- 
ered herein, no calculated T B was found. The behavior of these temperatures is 

"“max 

relevant to the second virial coefficient B and its relation to the principle of corre- 
sponding states (see appendixes A and B). 


Corresponding States Fluids 

In figure 8(a) the inversion curve for methane is in good agreement with the limited 
methane data of Tester (ref. 22), except near the knee of the curve. As the isenthalps 
are extremely flat in this region, small errors in PVT and in the partial derivatives can 
result in large shifts in the inversion curve. As the PVT data used here are from ref- 
erence 14, it is assumed that the inversion locus is more accurate than that of Tester 
(ref. 22). 

In figures 8(b) to (d) the nitrogen and oxygen inversion curves and data are replotted 
in reduced coordinates. The agreement has been discussed previously. 

In figure 8(e) the inversion curve for argon is plotted along with the data of Roebuck 
and Osterberg (refs. 4 and 7); Grossman, McCarty, and Hust (ref. 23); and Michels 
et al. (refs. 24 to 26); and the Kihara potential as calculated by Gunn (ref. 10). The 
agreement is quite good up to T B = 3. From T B = 3 to the maximum inversion tem- 
perature the relative error is at most 5 percent with respect to the results of refer- 
ence 23 . Since the calculated data points of reference 23 and the inversion curve of fig- 
ure 8(e) are both based on an extrapolation of the equation of state, the deviations seem 
quite reasonable. As such, the relative error could be used as an indicator of the con- 
fidence which can be placed in the maximum inversion temperature. 

In figure 8(f) the inversion curve for carbon dioxide is in good agreement with the 
data of Roebuck, Murrell, and Miller (ref. 27) and Newitt, Pai, Kuloor, and Huggill 
(ref. 28). Note that the shape of the reduced inversion curve deviates significantly from 
that of argon (fig. 8(e)). We also note that the acentric factor for argon is -0. 002 and 


o 

The quantum fluids, helium, hydrogen, and neon, are treated separately. Neon is 
considered a weak quantum fluid. 
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Figure 10. - Reduced Joule-Thomson inversion data for krypton from 
reference 44. Critical pressure, P c = 5. 502 MN/m 2 critical tem- 
perature, T c = 209.4 K. 
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Figure 11. - Reduced Joule-Thomson inversion data for xenon from 
reference 45. 



TABLE I. - FIXED- POINT TEMPERATURES RELATED TO 


SECOND VIRIAL COEFFICIENT B 


Fluid 

Maximum inversion 

Low- density Joule 

Boyle tempera- 

Reduced temperatures 


temperature 

inversion tern- 

ture (B = 0), 





(B/T =■ maximum), 

T. 

inv 

perature (B = 
maximum), 

t b 

max 

T B 

T R, inv 

TR ’ B max 

t r,b 

CH. 

1008.55 

— 

518.9 

5.2867 


2. 72 

V 

626.78 

2140 

321.4 

4.9627 

16.94 

2. 544 

°2 

798. 98 

— 

403.3 

5. 162 


2. 606 

a r 

815.06 

— 

412.7 

5. 4085 


2. 738 

co 2 

1387.78 

— 

725. 8 

4. 5619 


2.386 

Ne 

236. 03 

— 

122. 9 

5.316 


2. 767 

He 

55. 13 

196 

27.5 

10. 601 

37. 7 

5.292 

CO 

647.51 

— 

347.9 

4. 872 


2. 618 

b»r 

N 2 

644. 4 

2140 

329.8 

5.102 

16. 94 

2. 61 


a From ref. 14. 
^From ref. 13. 


TABLE II. - COMPARISON OF CRITICAL ISENTHALPIC JOULE- THOMSON COEFFICIENTS 


AS CALCULATED FROM REFERENCE 15 AND M c = T c /6P c 


Fluid 

Acentric 

Critical pressure, P c 

Critical tern- 

Critical 

Critical isenthalpic Joule-Thomson coefficient 


factor 
(ref. 50), 
w 

atm 

MN/m 2 

perature, 

T c 

density, 

Pc ’3 

g/cm' 3 

6P„ 

c 

K/MN/m 2 

LL ~- 
C 

(from 

/>P \ 

V'T / 

. (U 

ref. 

-1 

c 

15) 

ch 4 

+0. 013 

45. 66 

4. 627 

190.77 

0. 162 

6. 87 

t 

. 88 


“ n 2 

+ . 04 

33. 72 

3. 417 

126.3 

.3105 

6. 16 

c 

.91 


°2 

+ . 021 

50. 16 

5. 083 

154. 78 

. 4325 

5.08 

5 

. 18 


Ar 

-.002 

48. 01 

4. 865 

150. 7 

.531 

5. 16 

5 

.58 


co 2 

+ . 42 

72. 87 

7.3835 

304.21 

. 464 

6.87 

5 

.56 


Ne 

0 

26. 19 

2. 6537 

44. 4 

. 483 

2.79 

4 


CO 

+ . 041 

34.53 

3. 4986 

132.91 

. 2997 

6.33 

3 

.96 


He 

0 

2.245 

.227464 

5.2 

. 0693 

3. 81 

7.9 


s 

+ .04 

33.56 

3.4 

126. 2 

.3053 

6. 19 

6 

.22 


Xe 

+ . 002 

58 

5. 877 

289.75 

1. 1052 

8.22 

- 

— 


Kr 

-.002 

54.3 

5. 502 

209. 4 

.908 

6.34 

- 

— 


H 2 

0 

12. 77 

1. 294 

32. 984 

. 03077 

4.248 

- 

— - 



a From ref. 14. 
^From ref. 13. 
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that for carbon dioxide is 0. 42 (table n and appendix C), which means that carbon diox- 
ide deviates considerably from the principle of corresponding states. 

In figure 8(g) the inversion curve for carbon monoxide departs significantly from 
the data of Leah (ref. 29); however, the departure is not quite so significant for the data 
of Michels, Lunbeck, and Wolkers (ref. 30). The most troublesome aspect of these ex- 
trapolated results is the significant departure from the general form of the inversion 
curve for other corresponding states fluids. The results of Leah (ref. 29) and those of 
reference 15 were compared in depth, and are discussed in appendix E. There do seem 
to be systematic deviations in the results of Leah (ref. 29) and those of reference 15; 
however, the deviations presist into the regime where reference 15 is more accurate 
than reference 29. Thus, the inversion curve is questionable at high temperatures. 

The corresponding states form is definitely favored. 


Quantum Fluids 

Helium, hydrogen, and neon are quantum fluids and, as such, depart from a corre- 
sponding states treatment. However, neon appears to behave more like a normal fluid 
and is usually considered to be a weak quantum fluid. 

The inversion curve for helium and the data of Baehr (ref. 31), Keesom (ref. 32), 
and Hill and Lounasmaa (ref. 33) are plotted in figures 9(a) and (b). In general, the data 
are in good agreement with the curve; however, one data point of Baehr (ref. 31) seems 
to be out of line. The maximum inversion temperature of helium has been calculated by 
several researchers: Beattie (ref. 34) and Keyes (ref. 35), 45 K; Kirkwood and Keyes 
(ref. 36), 54 K; Roebuck (refs. 2 and 3) estimated 23. 6 K by using the reduced temper- 
ature of 4. 55 K for air; Gunn, Chueh, and Prausnitz (ref. 10), 47. 7 K; and an observed 
value of 34 K as listed by Crawford (ref. 37). Dean and Mann (ref. 19) use 55 K in their 
analysis of the Joule- Thomson refrigerator. The maximum inversion temperature of 
table H (appendix C) is 55 K. The agreement of the curve with reference 19 is not sur- 
prising as the equation of state used in reference 15 is that of Mann (ref. 38). The re- 
duced inversion curve for helium (fie:. 9(b)) is based on the pseudocritical parameters 
T c = 9. 4 K and P c = 0. 5492 MN/m as suggested by Gunn, Chueh, and Prausnitz 
(ref. 10). Superimposed on this figure is the helium inversion curve determined in ref- 

3 

erence 10 for comparison. Generally, the discrepancies are large on the upper part of 
the inversion curve, as would be expected from the lack of agreement in the maximum 


O 

Gunn, Chueh, and Prausnitz (ref. 10) developed three distinct curves: one for he- 
lium, one for hydrogen, and one for the corresponding states fluids. 
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inversion temperature, 47. 7 K (ref. 10) to 55 K (appendix A and table I). This would 
mean that the pseudocritical temperature would have to be 10. 8 Kin order to force 
agreement with the general reduced inversion temperature of 5. 095 K (ref. 10). 

In figure 9(c), the inversion curve for hydrogen was determined by using the pro- 
gram of Goldberg and Haferd (ref. 39), which is based on the PVT fit of Roder, Weber, 
and Goodwin (ref. 40). The data of Michels, DeGraaf, and Wolkers (ref. 41) are gener- 
ally in good agreement with the curve; however, they appear high. The lower tempera- 
ture data of reference 40 are in good agreement, as would be expected. In figure 9(d), 

the reduced hydrogen inversion curve is based on the pseudocritical values of refer- 

o 

ence 10, T c = 41.4 K and P c = 1.854 MN/m . As might be expected, the agree- 
ment between the curves of Gunn, Chueh, and Prausnitz (ref. 10) and the inversion 
curves is quite good as they came from the equivalent data sets. 

In figure 9(e) the reduced inversion curve for neon agrees very well with the 
smoothed data of McCarty and Stewart (ref. 42) and Michels and Gibson (ref. 43). This 
is not too surprising since the maximum inversion temperature was used as a constraint 
in determining the equation of state (ref. 42). In figure 9(f) the reduced neon curve is 
based on the pseudocritical values T = 45. 3 K and P = 2. 705 MN/m^ (ref. 10). 

In comparing the curves (figs. 9(e) and (f)) there appears to be little difference between 
them. (Recall that neon is considered a weak quantum fluid. ) 


Other Joule-Thomson Data 

The data of Trappeniers, Wassenaar, and Wolkers (ref. 44) and Michels, 
Wassenaar, DeGraaf, and Wolkers (ref. 45) were used to obtain portions of an inversion 
curve for krypton and xenon, respectively (figs. 10 and 11). The general trends of the 
data appear to be in good agreement with other fluids and are included in order to form 
a more general reduced inversion curve applicable to several fluids, as discussed later. 


Equations for Inversion Curves 

The reduced inversion curves for methane, nitrogen, oxygen, carbon dioxide, car- 
bon monoxide, argon, neon, helium, and hydrogen have been fitted by using a least- 
squares technique from Tr- 0.9 (P R ^ 2) to the maximum inversion temperature with 
a polynomial 

6 

P R =S A i T R 
i=0 
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TABLE III. - POLYNOMIAL COEFFICIENTS OF REDUCED INVERSION CURVES OF SEVERAL FLUIDS 


P R = 


E 

i=0 


A.T: 


i*R 


Fluid 

Number of 
points 
used in 
each fit, 
n 

A o 

A 1 

A 2 

A 3 

A 4 

A 5 

A 6 

Average sums 
of square of 
residuals, 
E r xl00 

Methane 

28 

-31.2651873 

60.3026109 

-32.9452939 

9.29996645 

-1.50083567 

0. 12645882 

-0.0041500129 

0.59 

a Nitrogen 

37 

-32.2267981 

64.2512980 

-37.7834487 

11.966361 

-2.29966244 

.25510747 

-.012830057 

.58 

^Nitrogen 

22 

-35.7879152 

72.8486137 

-45.6317248 

15.2554526 

-2.90682623 

.28560901 

-.010865648 

1.1 

Oxygen 

34 

-36.2893677 

75.0939713 

-48.741621 

17.2310336 

-3.57529548 

.40633931 

-.01973906 

.77 

Argon 

30 

-25.2735214 

47.0015602 

-22.2479188 

5.10306394 

-.67161214 

. 054807958 

-.0025724488 

.31 

Carbon 

26 

-49.6820235 

106.352075 

-72.1171474 

25.3427205 

-5.00402093 

.51606354 

-.021543775 

.9 

dioxide 










Neon 

30 

-41.3301625 

86.5178852 

-59.0821381 

21.5253196 

-4. 46121949 

.49273197 

-.02267271 

7.6 

Carbon 

32 

-24. 5141408 

44.1764731 

-18.5476859 

3.23409545 

-.26868672 

.014909515 

-.001050672 

1.5 

monoxide 










Hydrogen 

21 

-15.5988252 

26.0321395 

-9.7459013 

2.4207304 

-.49105816 

. 05932495 

.002913248 

. 18 

Helium 

26 

-9.50023901 

13.4661013 

-2. 1997754 

. 14814825 

-.01540754 

.00173746 

-.7093585x10' 4 

1 . 1 

General 0 


^0 

*1 

^2 

% 

^4 


^6 




-32.5209374 

65.6922312 

-39.7384300 

12.9300299 

-2.46176904 

0.25378553 

-0.010986565 



a From ref. 14. 

^From ref. 13. 

c For fluids with small acentric factors. Helium, hydrogen, and neon are not included due to quantum effects. 


The coefficients for these fluids are given in table HI, along with the average sums 
of the squares of the residuals 


E r n T"! ( P R ~ P R, data^ 

Generalized Reduced Inversion Curve 

In forming a generalized reduced inversion curve, inversion data for the fluids 
argon, methane, oxygen, and nitrogen and the data for fluorine (ref. 46), xenon, and 
krypton were weighted equivalent. The carbon monoxide data of Leah (ref. 29) and 
Michels, Lunbeck, and Wolkers (ref. 30) and points from the inversion locus of fig- 
ure 8(g) were weighted equivalent but not used in the fit; however, they are centered 
about the general inversion curve. The fluids helium, hydrogen, and neon were not used 
because of the quantum effects cited previously. Carbon dioxide was not used because 
of its large acentric factor (0. 42, table n and appendix C). Within these restraints, and 
recalling that much of the curves were calculated from extrapolated PVT data as dis- 
cussed previously, the locus of points in figure 12 is a least- squares representation of 
all the data for the corresponding states fluids. The generalized inversion curve is 
described by 


6 

P R=2>i T k 

i=0 

The coefficients /3. are tabulated in table m. The agreement between the general re- 
duced inversion curve presented herein and that of Gunn, Chueh, and Prausnitz (ref. 10) 
is quite good for T R < 3, as can be seen from figure 12. 

Although the quantum fluids and their isotopes should be handled separately (e. g. , 
refs. 10 and 21), the reduced inversion curves could probably be brought together by 
finding (empirically) the suitable pseudocritical temperature and pressure for each fluid. 
Figure 13 illustrates the pseudoreduced inversion curves for the quantum fluids relative 
to the generalized reduced inversion curve of figure 12. 
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o Methane (ref. 22) 

? Argon (refs. 4, 7, 10, 23 to 26) 

□ Oxygen (ref. 12) 

o Carbon monoxide (fig. 8(g)) 

o Nitrogen (refs. 5, 16 to 18, and 20) 

> Krypton (ref. 44) 

< Xenon (ref. 45) 

o Carbon monoxide (ref. 30) 

^ Carbon monoxide (ref. 29) 

^ Fluorine (ref. 46) 

Gunn's inversion curve (ref. 10) 

General inversion curve (this report) 

[ i T 1 1 T i 1 1 i 1 1 1 1 1 1 1 M 1 1 1 1 1 1 1 1 1 ii r 1 1 ii i n iim 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 m 1 1 1 1 1 in i 
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2 4 6 8 10 ] 

p R 

Figure 12. - Generalized reduced Joule-Thomson inversion curve for several corresponding states fluids, based on the PVT extensions 
and data. Curve of reference 10 included for comparison. 



Figure 13. - Generalized reduced Joule-Thomson inversion curve of 
figure 12 and pseudoreduced inversion curves of quantum fluids 
helium, hydrogen, and neon. Inversion curve for helium of ref- 
erence 10 included for comparison. (For nonquantum fluids, 

Tr " t R and p R = p R-» 


Joule-Thomson Coefficients 

With the aid of a computer program (e.g. , ref. 15), the isenthalpic Joule-Thomson 
coefficients (eq. (1)) can be calculated over a wide range of temperatures and pressures. 
Curves of p against T with P as a parameter, T against P with (jl as a param- 
eter, or various combinations can be fabricated. 


Reduced Joule-Thomson Coefficients 

In appendix C, a technique for determining the Joule-Thomson coefficient at the 
critical point is discussed 



(C4) 
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or, in normalized form, in terms of the dimensionless parameter, 


T c 

P c^c 



(C5) 


The difficulty is, of course, obtaining a proper derivative in the critical region; and 
based on the work of Bender (ref. 14), Coleman and Stewart (ref. 13), and reference 15, 
the normalized dimensionless parameter T c /P c p c is nearly a constant. Thus, a sim- 
plified form of p c is proposed 

T_ 

P r w — — (C7) 

C 6P c 


For the fluids discussed herein, the values of p c are tabulated (table II and appen- 
dix C). As can be seen the values are quite reliable except for helium, neon, and carbon 
monoxide. The equations of state for these fluids were fit many years ago and should 
probably be reexamined in the critical region. It must be pointed out that the values of 
P c for nitrogen as calculated by reference 15 from the state equations of Coleman and 
Stewart (ref. 13) and Bender (ref. 14) are 6.22 and 5.91, respectively. Moreover, 

Roebuck and Osterberg (ref. 5) experimentally determined (p ) = 5.72. The agree- 

V c/n 2 

ment with equation (C6) = 6. nj is quite good when the difficulties in calculating 

2 

are considered, 
c 

By using p = T /6P , reduced Joule -Thom son coefficients (p R = p/p ) as func- 
tions of reduced temperature T R with reduced pressure Ppj as the parameter are 
presented in figure 14 for methane, nitrogen (ref. 14), oxygen, argon, carbon dioxide, 
carbon monoxide, neon, and nitrogen (ref. 13). 


(3P/3T) p J 


Constant Joule-Thomson Loci 

Interpolating the p-T figure with P as the parameter, curves of constant p can 
be plotted as a function of temperature and pressure. The calculated values from the 
state equations of references 13 and 14 in reference 15 for nitrogen for constant p are 
compared with the data of Roebuck and Osterberg (ref. 5) in figure 15. As can be seen, 
the calculated values are in good agreement; however, there is some disagreement in 
the p < 0 regime between the calculated values and data of reference 5. 
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(a) Methane. 



(c) Nitrogen (based on ref. 131. 


Figure 14. - Reduced Joule-Thomson 






Figure 15. - Joule-Thomson coefficient contours based references 13 and 14, calculated from reference 15, as com- 
pared to the contours presented in reference 5. 


Joule-Thomson Refrigerators and Liquifiers 

Joule-Thomson refrigerators and liquifiers usually operate in the lowest tempera- 
ture portion of the cycle. The thermodynamics and operating optimums are discussed in 
appendix D. It was determined that the governing equations for the optimum figure of 
merit (compressor power input/refrigeration output) for both the refrigerator and the 
liquifier were of the same form. Further, if the amount of gas liquified were zero, the 
governing equations for the refrigerator and liquifier were identical. Dean and Mann 
(ref. 19) solved the governing equations for helium, hydrogen, and nitrogen. Their re- 
sults and those herein point out that maximum refrigeration (liquification) does not cor- 
respond to the maximum figure of merit, that is, optimum operating costs. 

A qualitative look at the Linde intermediate pressure cycle reveals that the optimi- 
zation procedure can only serve as a guide to the experiment. For example, the anal- 
ysis tends to indicate that the intermediate pressure should be near P c with rig/m > 
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0. 15. In the Linde process for liquifying air, the intermediate pressure is 40 atmos- 
pheres and n 2 /m = 0. 2 (ref. 47). 

The operating costs can also be reduced by precooling the inlet gas and/or cascad- 
ing using several single J-T cycles with different operating fluids. 


CONCLUSIONS 

The pressure- volume-temperature (PVT) relations for several fluids (methane, 
nitrogen, oxygen, argon, carbon dioxide, neon, carbon monoxide, fluorine, hydrogen, 
and helium) have been used to study the behavior of the isenthalpic and isothermal Joule- 
Thomson coefficients and in particular the inversion curve. The following conclusions 
were reached: 

1. Examination of the PVT data obtained by extrapolating the equation of state 
described by Hendricks, Baron, Peller, and Pew indicates valid results for the high- 
pressure and high -temperature ranges. This is further validated by the fact that the 
partials (0P/9p),p and (3P/9T) , where p is density, are in good agreement with other 
high-pressure and -temperature results of Coleman and Stewart. The Joule -Thom son 
inversion curves calculated as functions of these extrapolated PVT data agree well with 
best available J-T data. 

2. Since the Joule- Thomson inversion curve (isenthalpic J-T coefficient p = 0) 
conforms to the principle of corresponding states, a generalized, reduced inversion 
curve can be determined for the fluids with low acentric factors. Although the individual 
fit of the inversion curve for each fluid gives better results, the generalized curve is 
useful in the absence of data and to indicate similarities among the fluids. The inver- 
sion curve developed herein and that presented by Gunn, Chueh, and Prausnitz are in 
good agreement. The quantum fluids helium, hydrogen, and neon are best handled 
separately. 

3. A critical isenthalpic Joule- Thomson coefficient p £ can be determined and a 
simplified expression deduced for it: 



Even though the corresponding states principle does not apply to p * 0, reduced curves 
can be determined and are useful in grouping data and where data are not available. 

4. The maximum inversion temperature (maximum (B/T)) and the Boyle tempera- 
ture (B = 0) can be found from the second virial coefficient B of the state equations. 
For the fluids with inversion curves calculated herein, no maximum in B was found. 


33 



5. Joule- Thomson liquifiers and refrigerators operate most efficiently at pressures 
substantially less than those corresponding to the maximum cooling capacity. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 10, 1972, 

113-31. 
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APPENDIX A 


PRINCIPLE OF CORRESPONDING STATES AND THE INVERSION CURVE 


The corresponding states principle applies to molecules with spherical symmetry 
and a universal curve shape for the potential energy of molecular interaction. Conse- 
quently, the principle applies only to those properties which depend on the intermolec- 
ular energy. As such, the reduced second virial coefficient B p Q should be a universal 
function of reduced temperature T R . Simple substances such as the inert gases (ex- 
cluding helium, hydrogen, and possibly neon due to quantum effects), nitrogen, oxygen, 
carbon monoxide, and methane do conform as expected. 

The second virial coefficient, to which the principle applies, is simply related to 
the isothermal Joule- Thomson coefficient 


in the limit of zero pressure 





lim 

P-0 



(Al) 


(A2) 


Consequently, the principle of corresponding states is applicable to cp . Now <p is 
related to the isenthalpic Joule- Thomson coefficient ju by 


<P = 




(A3) 


Francis and Luckhurst (ref. 48) note that the principle of corresponding states in 
the low-pressure limit is not applicable to specific heat C or to the isenthalpic Joule- 
Thomson coefficient ji, even though their product is a corresponding states variable. 

However, the application of the principle of corresponding states to the inversion 
curve, for the fluids indicated previously, is valid. The inversion curve represents the 
locus of points where 



(A4) 
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This, of course, also means that isothermal Joule- Thomson expansions for <p < 0 
produce cooling and isothermal expansions where cp > 0 produce heating. However, 
the inversion curve, common to both isothermal and isenthalpic expansions, represents 
a locus of minimums of the isotherms in the enthalpy- pres sure plane (fig. 16). 


Temperature, 

K 
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APPENDIX B 


INVERSION TEMPERATURES AND THE SECOND VIRIAL COEFFICIENT 


The isothermal and isenthalpic Joule- Thomson coefficients, as discussed in appen- 
dix A, have been used to test virial forms of the equation of state and derived proper- 
ties. 

From reference 15 the equation of state can be written as 
P = pRT + A(T)p 2 + B(T)p 3 + C(T)p 4 + D(T)p 5 + H(T)p 6 


[e(T) + F(T)p 2 + G(T)p 4 ] 


3 ( -n 20^ Z^o) 

p 6 e Z (Bl) 


where the functions H(T) and G(T) were added to accommodate Coleman and Stewart's 
equation for nitrogen (ref. 13). Expanding the exponential in a Taylor series about 
p = 0 and rearranging, the equation of state can be written as 

£ = z = 1 + Ail! P + B(T) + E(T) e ( “ n 20 p2/T 2 ) p 2 + C(T) p 3 + _ _ 
pRT RT RT RT 


= 1 + Bp + Cp 2 + Dp 3 + . . . (B2) 

where the first term represents a fluid with negligible molecular interaction (i.e. , ideal 
gas). The second B, third C, or fourth D, etc., terms or virial coefficients repre- 
sent interactions between two, three, or four molecules, etc., respectively. As such 
the virials B, C, D, etc. , represent the deviation from an ideal gas state. Here we 
are interested in the second virial coefficient B = A(T)/RT, or from reference 15 with 
the modifications cited. 



(B3) 


The zero- pressure inversion temperature (maximum inversion temperature) is 
given in terms of the second virial coefficient for the equation of state 


PV = RT 1 +— + — + . 


= 1 + B'P + C'P^ + 


(B4) 


37 



by the zero-pressure limit of <p (appendix A). Using equations (Bl) and (1) we find the 
limit 


Urn (-< -p) = lim liC = T 2 (B5) 

P-0 P-0 p dT 

Since C p > 0, ju passes through zero at the temperature for which B/T is a max- 
imum. By using the second virial coefficient of reference 15 (eq. (B3)) and maximizing 
B/T, the zero-pressure inversion temperatures T iny estimated by extrapolating the 
PVT results to higher temperatures are presented in table I. For the Van der Waal 
equation of state 


P RT _ _a_ 
" V - b y 2 


(B6) 


the second virial coefficient becomes 


B = b — — (B7) 

RT 

and substituting equation (B7) into equation (B5) the maximum inversion temperature 
becomes 


T. 

mv 


2a 

Rb 


(B8) 


or in reduced form 


T 

R, inv 


27 

4 


(B9) 


since Van der Waal's value of T c is 8a/27Rb. Meissner in the Handbueh der Physik 
(1926) found T R ^ ny to be about 5. 8 for all gases, which is 14 percent lower than given 
by equation (B9) except, of course, for helium. However, the values presented in 
table I are even lower but are generally in good agreement with the results of Gunn, 
Chueh, and Prausnitz (ref. 10), where T R ^ ny is about 5.075 K. For nitrogen, ref- 
erence 34 gives T r ^ = 4. 75 and reference 6 gives T R iny = 4. 92. For argon, ref- 
erence 6 gives T r iny = 4.8, and reference 23 gives T R - ny = 5.26. 

The Joule or free expansion coefficient 
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V s P 


2 



has an inversion temperature at the maximum in the second virial coefficient B for 
low densities. While the relation between the Joule-Thomson coefficients <p and /j. 
is simple, the relation to the Joule coefficient is more complex 

’ c * =T (i) p - p 

= P 2 1- 9 * 2 

where <p = - pC . Substituting equation (B4) into (BIO), we find for the region of low 

r 

pressure (low density) that 



77 C v _ dB 
Rp 2 T 2 dT 


(Bll) 


Note that at zero pressure, 77 is zero; but for densities low enough for equation (Bll) 
to hold, the Joule coefficient has an inversion where B is a maximum. Values of this 
inversion temperature, as listed in table I, were obtained from reference 34. B as 
given by equation (B3) does not possess a maximum for T ^ 4000 K. Such tempera- 
tures are far beyond the range of the curve fits and beyond the scope of the Joule- 
Thomson curves for fluids discussed herein. It is also interesting to note that this in- 
version temperature does not exist for the Van der Waal fluid, as B does not attain a 
maximum . 

The Boyle temperature is that temperature for which the second virial coefficient 
vanishes, that is, when B = 0. The Boyle temperatures are also given in table I. For 
the Van der Waal fluid, equations (B8) and (B9) give 


T R, B 


T 

R, inv 
2 


(B12) 


The relation between calculated and reduced temperatures and the calculated reduced 
Boyle temperatures seem to follow equation (B12) even though the Van der Waal tem- 
peratures Tg and T iny do not. 
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APPENDIX C 


CRITICAL ISENTHALPIC JOULE -THOM SON COEFFICIENT 

The specific heat at constant pressure may be written as 



When this expression for C is substituted into equation (1), the isenthalpic Joule- 

r 

Thomson coefficient becomes 



Now as one approaches the critical point, the equation of state strongly satisfies 

lim (^\ - 0 (C3) 

PVT-P c v c T c Va P ; T 


while the specific heat at constant volume has at most a logarithmic infinity; that is, 
very slowly approaches 00 . Thus, the critical isenthalpic Joule- Thomson coefficient 
becomes 


M 


c 


lim 

PVT ~P C V C T C 


M = 



PVT 
c c c 


(C4) 


Furthermore, multiplying equation (C2) by P/T, using the result of equation (C3), 
and applying the limit of equation (C4) yields the dimensionless parameter 
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(C5) 


P c^c _ 

lim PH = 

/a In p\ 

-1 

T c " 

PVT-P C V C T C T 

U In Tj p _ 



which normalizes the critical isenthalpic Joule-Thomson coefficient p c . The value of 

(0P/0T) p can be found from references 14 and 15. However, for methane, nitrogen, 

oxygen, carbon dioxide, and argon Bender (ref. 14) gives a nearly constant value for the 

product (T/P)(0P/0T) p v T ; that is, the normalized critical isenthalpic Joule- 

*c v c 1 c 

Thomson coefficient (eq. (C4)) is nearly a constant, 


T c 

P c^c 


- “ 


- 

/a In p\ 

_ ^C 

(*£' ) 

\0 In T I 

P_ 

\ 9T L 

\ f Q 

PVT 

c c c 

x / p 


(C6) 


Rearranging \± c becomes 


n k—Z. (C7) 

c 6P 

c 

The values of ji Q calculated from equation (C7) and from reference 15 for the fluids 
considered herein are given in table II . 

Using i± c we can find the critical dependency of the isothermal Joule-Thomson 
coefficient cp and the Joule or free expansion coefficient 77. The Joule or free expan- 
sion coefficient may be written as 


^ C v 



- P 


Rearranging and passing the limit 


lim I 

PVT— P c V c T c 


+ PN 


lim 

pvt-p c v c t c 


t/0P\ 
P \ 3T / 
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and 





since ^C y j is weakly divergent, rj c must approach zero in a weak manner. 

The isothermal Joule- Thomson coefficient, from appendix A, is directly related to 
li through C , 

r 


(p = - MC 


In the critical region, C diverges as |t - T c |" 4//3 , while 11 remains finite and non- 


zero. Consequently, cp diverges as C , that is, 


<P. 


j T - T j ” 4 ^ 3 
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APPENDIX D 


JOULE-THOMSON LIQUIFIER AND REFRIGERATOR 
CYCLES AND OPERATING OPTIMUMS 

A Joule- Thomson liquifier and refrigerator is illustrated schematically in figures 17 
and 18 (see also refs. 19 and 49). The only difference between the two devices is the 
makeup gas supplied to the liquifier at the compressor; consequently, the temperature- 
entropy and enthalpy-temperature diagrams apply to either process. 



Figure 17. - Process diagram and flow schematic for a Joule-Thomson liquifier. Process diagram 
also applies to Joule-Thomson refrigerator. 


The compressor need not be a part of the cycle; however, high-pressure gas must 
be supplied and a return provided. The Joule-Thomson process could be the last of sev- 
eral stages in a real refrigerator or liquifier of gas compressed isothermally at ambient 
temperature T a to inlet pressure P^. 

Upon passing through the heat exchanger m g/sec of fluid is expanded isenthalpi- 
cally such that n g/sec of the gas are liquified. The energy balance requires 

riiHj = nH 4a + (m - n)Hg (Dl) 


Since we want to maximize the amount of liquid produced and the liquid reservoir 
pressure is fixed, the only variable at our disposal is H^, which is controlled by P-^ 
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Figure 18. - Process diagram and flow schematic for Joule- 
Thomson refrigerator. Process diagram also applies to Joule- 
Thomson liquifier. 



(D2) 


But the latter term is the isothermal Joule- Thomson coefficient. This term is zero 
only when P and T lie on the inversion curve (appendix A). Consequently, the maxi- 
mum liquification or refrigeration occurs when the process starts on the inversion 
curve . 

However, the maximum liquification or refrigeration "costs" more than operating 
on a locus within the Joule- Thomson inversion curve. To see this, again consider the 
liquifier. If we let rj e represent the efficiency of the compressor, etc. , the work re- 
quired to compress m g/sec of gas (isotherm ally) becomes 

W„ = riiRT In 

P 3- 

Dividing by the amount of liquid produced (n = fig g/sec) gives a figure of merit for the 
system. Minimizing the figure of merit with respect to our only variable gives 
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and from the heat exchanger load, 

(Hj - H 2 ) = ^1 - ^ (Hg - H 4f ) (D4) 

which represents the second law constraint (heat flow requires a positive temperature 
difference). Also, it must be noted that the process is constrained to those points where 
Pg is less than P c and q is less than zero. 

Equation (D3) is identical in form to that of Dean and Mann (ref. 19) for Joule- 
Thomson refrigerators. Note that equation (D4) for the Joule- Thomson liquifier readily 
reduces to the heat exchanger load for the refrigeration, n 2 equals zero. Dean and 
Mann numerically solved equation (D3) with n 2 replaced by m and equation (D4) with 
[l - (n 2 /m)] replaced by unity, for the fluids helium, para- hydrogen, and nitrogen, pre- 
senting curves similar to that of figure 19. It is quite apparent that one should operate 
well within the inversion curve to minimize cost. However, other methods of reducing 
the operating costs are to add intermediate cycles, as illustrated in figure 20, to pre- 
cool the inlet gas, or by cascading using second single cycles with different fluids. 



Inlet pressure, Pj 

Figure 19. - Comparison of inversion curve and 
optimum (minimum costloperating locus for 
Joule-Thomson processes. 
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1 3 5 



Figure 20. - Double-cascade Joule-Thomson liquifier 
(early Linde process for air). 


The energy balance, isothermal compressor work, and heat exchanger balance be- 
come, for ju > 0, 


mH 1 = n 2 H 2f + (m - n 2 )H 3 
n 2 H 2f = n 4 H 4f - (n 2 - n 4 )H 5 



m(H 1 - H a ) = (m - n 2 )(H 3 - Hg) + (n 2 - n 4 )(H 5 - H 4 ) 


Note here that (m - n 2 ) represents the mass of gas returning to the compressor and that 
n 2 represents the gas or liquid passing to the second Joule-Thomson valve. Solving 
for the figure of merit, 


46 




The optimization procedure becomes quite complicated and of academic interest. 

For example, the optimum Wp/m in equation (D6) is for rig equals zero; that is, no 
fluid passing to the second J-T valve. However, the whole idea is to produce liquid, and 
the optimum curves are found by plotting those points that can be attained by expanding 
from Pj to Pg. The analysis can, however, serve as a guide to experimental testing. 

If one reflects on the previous discussion for a moment, and uses a little intuition 
along with figure 7, one would anticipate the maximum temperature drop to occur if the 
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process started at the inversion curve. Furthermore, since most of the curves are 
quite flat near the inversion curve, one would expect that nearly the same amount of 
refrigeration could be obtained for a significantly lower pressure drop, at a substantial 
savings in compressor power. 
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APPENDIX E 


DIFFERENCES IN CARBON MONOXIDE PROPERTIES 

The inversion curve of figure 8(g) as calculated from reference 15 and the results of 
Leah (ref. 29) are in considerable disagreement. As the PVT data from reference 15 
are extrapolated considerably beyond the range of data used in the curve fit (Tr < 2^, 

Pr < 9), a direct comparison of density, enthalpy, entropy, and specific heats was made. 

In figure 21(a) the difference in density is plotted as a function of reduced pressure 
for temperature increments of 100 K. For T R < 3, the difference appears to be about 
0. 2±0. 1 percent, which seems quite good. The shaded band indicates the difference ex- 
pected at the inversion curve for Tr > 2. One part of the band is formed from fig- 
ure 8(g), the other from the data of Leah (ref. 29). As can be seen, the difference in- 
creases as Tr increases to 5 and Pr ^ 6. The difference in density along the 
inversion curve then diminishes as pressure is reduced. Thus, the inversion curve in 
the region Tr > 3 and 1 < Pr <11^ seems to have a significant departure from that 
of Leah (ref. 29). However, in figure 8(g), the deviations from Leah (ref. 29) are sig- 
nificant even for Tr < 3. Thus, not all the answers lie in the PVT results. 

In figure 21(b) the difference in enthalpy along the inversion curve is considerable 
between 10 < Pr < 12, passes through a minimum near Pr = 9, and increases as pres- 
sure decreases. As very small deviations in the enthalpy in this region will significantly 
alter the inversion locus, it appears that derived properties (i. e. , the (3P/3T) ) may be 
the source of the differences. 

However, in figure 21(c) the difference in entropy along the inversion curve is about 
0.9±0.2 percent for Pr < 11, which apart from the level shift in Sq is quite reason- 
able. As both enthalpy and entropy require (3P/3T) , such agreement tends to exoner- 
ate the derivative (3P/3T)^. 

In figure 21(d) the difference in specific heat at constant volume is quite large and 
directly related to temperature level. Even in the regime where the results of refer- 
ence 15 are quite accurate (Tr < 2^) the differences are large. This, however, would 

2 2 

only tend to implicate 3 P/3T , which is a troublesome point in many a curve fit and 
sheds little light here, as either (or both) references 15 or 29 could be in error. 

In figure 21(e) the differences in specific heat at constant pressure along the inver- 
sion curve are, as were C y errors, a function of temperature level. However, super- 
position of the two difference curves indicates them to be equivalent, except at Tr < 3 . 
As C is determined from C and the partials (3P/3T) and (SP/SpL, seem valid, 

p v 2 7 2 1 

C appears to be within the limits of uncertainty in (3 P/3T ). 

P 

Consequently, it would appear that the variations in density coupled with the varia- 
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Reduced 

temperature, 

T R 


Temperatu re, 
T, 

K 


° 2.055 273 

a 2. 806 373 

o 3. 56 473 

* 4.31 573 

-a 5.065 673 


CO inversion from data (fig. 8(g)) 

(ref. 29) 


CO inversion locus (fig. 8(g)) 



(c) Relative entropy differences. 

Figure 21. - Comparison of carbon monoxide properties as calculated from reference 15 with those of reference 29, at 
several reduced pressures for selected isotherms. Expected differences along inversion curves (Leah's results vs. 
ref. 15) are bounded by shaded region. 
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tions in (0P/9T)^ lead to small shifts in the isenthalps, which in turn lead to deviations 
in the inversion curve. However, the data of Michels, Lunbeck, and Wolkers (ref. 30) 
depart significantly from those of Leah (ref. 29) and to a lesser extent from the results 
of reference 15. But these data appear to follow a curve about 5 percent higher than the 
inversion curve of figure 8(g) and 10 percent lower than that of Leah (ref. 29). 

It is quite apparent that the carbon monoxide data need to be reevaluated and a new 
equation of state presented. 

Finally, invoking the principle of corresponding states, one would expect the inver- 
sion curve to more closely follow that of Leah (ref. 29). Therefore, the curve of fig- 
ure 8(g) should be considered valid for < 3 and low for > 3. 
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APPENDIX F 


SYMBOLS 


A(T), B(T), C(T), D(T), 
E(T), F(T), G(T), H(T) 

A 0,l,. • -6 
a,b 

B 

B* 

C 

C' 

S 

C v 

D 

H 

k 

m 

N 


P 

R 

S 0 

U 

V 

w 

p 

T 

T 2 


temperature functions in equation of state of ref. 15 
constants 

constants in Van der Waal equation of state 

second virial coefficient 

modified second virial coefficient 

third virial coefficient 

modified third virial coefficient 

specific heat at constant pressure, J/(g)(K) 

specific heat at constant volume, J/(g)(K) 

fourth virial coefficient 

enthalpy, J/g 

summation index 

mass flow rate, g/sec 

number of stages 

number of points 

mass flow rate, g/sec 

constants in equation of state (ref. 15) as modified to accommo- 
date ref. 13 (i = 1, 2, 3, 4, 5 . . . 20, 21, 25) 

9 

pressure, MN/m 
gas constant, J/(g)(K) 
reference entropy, J/(g)(K) 
internal energy, J/g 

O 

specific volume, cm /g 
compressor or pump work, J/(sec)(W) 
temperature, K 

temperature, used in helium fit only, K 
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z 

a i 

*1 

p 

<p 

n 

% 

p 

00 

Subscripts : 
a 
B 
c 

data 

inv 

f 

max 

N 2 

R 

vp 

0 

1,2,3 . . . 


compressibility factor, PV/RT 
constants (i = 1, 2, ... 5) 
constants (i = 1, 2, . . .7) 

q 

density, g/cm 

q 

isothermal Joule- Thomson coefficient, cm /g 

q 

Joule or free expansion coefficient, (g)(K)/cm 
efficiency 

isenthalpic Joule- Thomson coefficient, K/(MN)(m ) 


acentric factor, 




(ref. 50) 


ambient 

Boyle 

critical 

data 

inversion 

liquid 

maximum 

nitrogen 

reduced 

vapor pressure 

zero pressure 

states on thermodynamic diagram 


Superscripts: 


0 pseudovalue 

’ (see fig. 9) 
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